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By Gaylord 0. Ellis and John D. Stanitz 

SUMMARY 

A solution is presented for three-dimensional, incompressible, non- 
viscous, potential flow in a rotating impeller passage with zero through 
flow. The solution is obtained for a conventional impeller with 
straight blades but with the inducer vanes removed and the impeller 
blades extended upstream parallel to the axis of the impeller. By super- 
position of solutions two additional examples are obtained for different 
ratios of compressor flow rate to impeller tip speed. The three- 
dimensional solutions are compared with corresponding two-dimensional 
solutions and it is concluded that, at least for the type of impeller 
geometry investigated, two-dimensional solutions can be combined to 
describe the three-dimensional flow in rotating impellers with suffi- 
cient accuracy for engineering analyses- 


INTRODUCTION 

As an aid to better understanding of flow conditions in rotating 
impeller passages, methods of analysis have been developed in the past 
for potential nonviscous flow. In order to achieve solutions with a 
reasonable expenditure of effort, all methods are based on two- 
dimensional assumptions, in that the flow is restricted, by assumption, 
to specified flow surfaces in space. Either of two types of surface 
are usually assumed for the flow: first, the mean blade (or passage) 

surface on which flow conditions vary from hub to shroud but are con- 
sidered constant in the circumferential direction (axial- symmetry solu- 
tions, references 1 and 2), or, second, surfaces of revolution on which 
flow conditions vary from one blade to the next, but normal to which 
the flow conditions are considered constant (blade -to -blade solutions, 
references 3 and 4). 

If the streamlines of an axial symmetry solution are used to 
generate surfaces of revolution around the axis of the impeller, the 
totality of the blade -to -blade solutions on these surfaces of revolution 
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constitute a quasi -three -dimensional solution (reference 5) "because the 
solutions indicate variations in flow conditions throughout the impeller 
passage. However, because the flow is constrained to surfaces of revo- 
lution, the solution is not three dimensional in the exact sense of the 
word. No complete three-dimensional solutions for rotating impeller 
passages exist in the literature, and a solution has therefore been 
obtained at the NACA Lewis laboratory. The solution is presented in 
this report and is compared with the results of axial-symmetry and 
blade -to -blade solutions in order to evaluate these two-dimensional 
methods of analysis. 

The three-dimensional solution was obtained for incompressible non- 
viscous flow in a rotating impeller passage with straight blades and 
with the inducer vanes located far upstream of the impeller. By super- 
position of solutions, results are obtained for several ratios of flow 
rate to impeller tip speed. 

GENERAL METHOD OF ANALYSIS 

A partial differential equation for three-dimensional flow in a 
rotating impeller passage is developed from considerations of continuity 
and absolute irrotational fluid motion. 

Assumptions. - The fluid is assumed to be inviscid and incompres- 
sible"! The flow is assumed to be steady relative to the rotating 
impeller passage, and in the absence of viscosity the absolute motion 
of the fluid is assumed to be irrotational. It is assumed that the 
phenomenon being investigated, that is, the deviation of three- 
dimensional flow from the restricted motion of two-dimensional solu- 
tions, is qualitatively the same for compressible and incompressible 
solutions. This deviation is a perturbation resulting primarily from 
rotation of the impeller ; and in reference 6 it is shown that at least 
for two-dimensional solutions this type of perturbation is independent 
of compressibility, which affects only the average velocity. 

Cylindrical coordinate system and velocity components . - The cylin- 
drical coordinates R, 0, and Z relative to the impeller are shown in 
figure l(a). (All symbols are defined in appendix A.) These coordi- 
nates are dimensionless, the linear coordinates R and Z having been 
divided by the impeller tip radius (so that R is equal to 1.0 at the 
impeller tip) . 

The absolute- velocity Q has components Qr, Qg, and i n "the 

R, 0, and Z directions, respectively (fig. 1(a)). These velocities 
are dimensionless, having been divided by the impeller tip speed (so 
that, for example, the dimensionless blade speed at any radius R is 
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equal to R) . If W Is the velocity of the fluid relative to the impeller, 
expressed as a ratio of the impeller tip speed, then 

W 0 = Q 0 - R (1) 

Potential function qp . - For absolute irrotational fluid motion 

VXQ = 0 (2) 


where the bar indicates a vector quantity. A potential function cp 

satisfies equation (2) identically if defined by 

Q = Vcp ( 3 ) 

from which 

|r = % ( 3a ) 


and 




(3b) 


(3c) 


Differential equation of flow . - From continuity 

V.Q = 0 


so that, from equation (3), 


V 2 cp = 0 


which in cylindrical coordinates becomes 


(4a) 


d 2 cp 1 dcp 1 _ afcp Sfcp 

dR 2 R r2 ^ q 2 ^ z 2 


(4b) 


NUMERICAL PROCEDURE 

A numerical procedure is outlined for the solution of the partial 
differential equation (4) for flow in a rotating impeller passage with 
special type of geometry. 
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Preliminary Considerations 

Special type of impeller geometry. - The three-dimensional solutions 
presented are for a straight-bladed impeller of conventional design 
except that the inducer vanes are removed and the straight impeller 
blades are extended indefinitely upstream parallel to the axis of the 
compressor. This idealized entrance condition along with straight 
blades results in substantial simplification of the numerical procedure. 

Superposition of solutions. - As a result of the special type of 
impeller geometry just discussed, the boundary conditions for flow 
through a rotating impeller axe equal to the sum of the boundary condi- 
tions for zero flow through the rotating inpeller and for finite flow 
through the stationary impeller. Therefore, because the boundary con- 
ditions can be added and because the differential equation (4) is 
linear, the velocity potential cp for flow through the rotating impeller 
passage can be expressed as 

Cp = cp. + kcp 2 (5) 

where cpp satisfies equation (4) and the boundary conditions for the 
rotating impeller with zero net through flow and satisfies the same 

equation but for the boundary conditions associated with flow through 
the stationary impeller. The solution for cp-^ is called the "eddy-flow 

solution" and corresponds to ideal flow conditions in the rotating 
impeller with the throttle closed so that no through flow occurs. The 
solution for cpg is called the "through-flow solution" and, for the 

special type of Impeller geometry being considered, this solution is 
axially symmetric and corresponds to flow with zero whirl through an 
annulus with the same hub-shroud profile and no impeller blades. Solu- 
tions for various ratios of flow rate to impeller tip speed are obtained 
directly for various values of k in equation (5). 


Eddy-Flow Solution 

The eddy-flow solution for the rotating impeller passage with zero 
net through flow is considered first. 

Transformation of coordinates . - It is convenient for purposes of 
the numerical solution by relaxation methods to transform the RZ-plane 
to one on which the coordinates are the strea m lines t] and velocity 
potential lines £ for flow through the compressor annulus without 
blades. Because the hub and shroud contours are streamlines in the 
RZ-plane, these contours become straight parallel lines in the 
1 1 ] -plane. In terms of the new transformed coordinates, equation (4b) 
for the eddy-flow potential cp^^ becomes (appendix B) 
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Qc 


a 2 ^ 
a | 2 


i a 2 cp L ]_ a 2 ^ 


aq>L 


+ r2Q 2 2 ~ + ^2 ^2” + 2(Q 2 } Z w = 


a n 


( 6 ) 


where the subscript 2 refers to the solution for axially symmetric 
flow through the compressor annulus with no hlad.es or , which is the 
same thing, through the stationary impeller passage of the special type 
considered in this report. 

The new coordinate system introduces two additional velocity compo- 
nent s ( appendix B ) 




dqp 

= Q2 

(7a) 


and 

% = RQ 2 ^ 

(7b) 


and an angle a ,2 

defined by 


- 


tana 2'Qf) 2 

(8) 

-ft 

all of which are 
that 

shown in figure l(b). From this figure it is 

seen 



= Qg cos 02 + Or sin 

(9a) 


and 

= Q r cos 02-0% sin ^ 

(9*>) 


or, conversely. 

Qr = sin ag + cos ag 

(10a) 


and 

0% = cos ct 2 - sin 02 

(10b) 


Boundary conditions . - For the eddy-flow solution of equation (6) 
the boundary conditions that must be satisfied for the special type of 
impeller geometry considered in this report are : 
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(l) The flow direction must he tangent to the huh and shroud in the 
impeller and diffuser so that is zero, or, from equation (Tb), 



(2) Along the blade the relative flow is tangent to the blade sur- 
face so that for straight radial blades the relative tangential velocity 
Wo is zero and from equations (l) and (3b) 

= R 2 

(3) Boundaries are established in the diffuser on meridional planes 
extending from the blade tips. For a rotating impeller with no through 
flow the radial velocity component is zero on these boundaries so that 
the potential function is constant along radial lines on these surfaces. 
Variations in velocity potential cp with Z at the impeller tip of 
constant radius indicate the presence of a vortex sheet shedding from the 
trailing edge of the blade and passing downstream. It is assumed that 
the strength of this sheet is weak and can be ignored in the solution of 
equation ( 6) . For impeller blades with constant tip radius the variation 
in work input from hub to shroud at the impeller tip is negligible and 
the assumption therefore appears to be reasonable. The Joukowski condi- 
tion at the blade tip is automatically satisfied by condition (2). 

(4) The domain of the solutions is extended in the upstream and 

downstream directions until flow conditions are uniform in a plane normal 
to the direction of through flow. For the eddy-flow solution this con- 
dition is achieved when (Qi)| is zero, that is, when is zero, 

everywhere on a plane normal to the £ coordinate. 

(5) The idealized inlet of the special impeller geometry considered 
in this report results in symmetry of flow about the mean plane between 
blades in the rotating impeller with no through flow. The flow is 
directed normal to this plane and CfXj is therefore everywhere constant 

(zero) on it. 

Relaxation solution . - The differential equation (6) is solved by 
relaxation methods (reference 7) to satisfy the boundary conditions just 
described. The velocity components are then determined by equations (3) 
and (7) in finite difference form. For the numerical examples of this 
report, a three-point system was used for expressing the differential 
equations in finite difference form. 
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Flow paths . - Any three velocity components determine the flow direc- 
tion at a point so that flow paths relative to the impeller passage can 
he determined from the velocity components Q R , Wg, and Qg, or Q^, W q , 

and Qyj . On the hub, shroud, and blade surfaces the path lines can be 
constructed graphically from lines of constant flow direction on these 
surfaces. 

Accuracy . - For the numerical examples of this report, the impeller 
channel includes a total of 5400 grid points at which the velocity poten- 
tial was relaxed to a unit change in the fifth decimal. (Because cp^ 

is constant on the mean plane and the flow is symmetrical about this 
plane, the number of grid points at which it was necessary to relax is 
reduced to 2400. ) 

In order to check the accuracy of the graphical construction of the 
path lines, these lines were obtained on a plane normal to the through- 
flow direction far upstream of the impeller where a direct two-dimensional 
solution for the stream function is known and valid. Figure 2(a) compares 
the path lines with the streamlines. It is noted that the graphically 
constructed path lines agree well with the streamlines. It should be 
pointed out, however, that the path-line spacing is not sufficiently 
accurate to be indicative of the velocity distribution. In figure 2(b) 
the velocities obtained from the three-dimensional solution for the 
velocity potential are compared with the velocities obtained from the 
two-dimensional solution for the stream function. The comparison indi- 
cates much better agreement in the velocity distributions than was indi- 
cated by the path-line spacing in figure 2(a). 

A check on the accuracy of the three-dimensional solution will be 
given in connection with a discussion of the numerical examples. This 
check indicates approximately the same accuracy that is shown by the 
comparison of velocities in figure 2(b). 


Combined Solutions 

After the eddy-flow solution has been obtained, various percentages 
of a through-flow solution may be added to obtain solutions for different 
ratios of compressor flow rate to impeller tip speed. 

Through-flow solution . - The through-flow solution is obtained by 
methods outlined in reference 1, for example. As already discussed, the 
velocity potentials for the two types of solution can be added or, as 
indicated by partial derivatives of cp in equation (5), the velocity 
components themselves can be added directly. The latter procedure avoids 
the necessity of computing the distribution of from the distribution 

of stream function determined by reference 1. 
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Flow path . - The procedure for graphically determining the flow path 
for the combined solutions is identical with that outlined for the eddy- 
flow solution. 


NUMERICAL EXAMPLES 

Three-dimensional solutions for flow through an impeller with 
straight blades and with the inducer vanes located far upstream of the 
impeller are presented for: (l) zero flow through the rotating impeller 
passage, (2) flow through the stationary impeller passage, and (3) com- 
binations of (l) and (2) for various ratios of through flow to impeller 
tip speed. 


Impeller Geometry 

The impeller geometry for the numerical examples is the same as 
that in references 1 and 4 with the inducer vanes located far upstream 
of the impeller. The hub-shroud profile of the impeller is described 
in figure 3. The blade spacing is 32.80° as in reference 4. 

The results of the solutions are presented on the channel surfaces 
and on the nine meridional planes indicated in figure 4(a). The |,T] 
coordinates on the meridional planes are shown in figure 4(b) . The lines 
of constant | are spaced at intervals corresponding to equal increments 
of the | coordinate used in reference 4. 


Solution for Zero Net Flow Through Rotating Impeller Passage 

Velocity potential cp-j^. - Lines of constant velocity potential on 

the meridional planes are shown in figure 5. The center plane E (see 
fig. 4(a)) is not shown because, as discussed previously, cp x is zero 

everywhere on this plane. Note that lines of constant cp-^ intersect the 

hub-shroud profile at right angles, as required by equation (ll). The 
meridional velocity component must be directed normal to the lines of 
constant in the meridional planes and has magnitudes inversely pro- 

portional to the line spacings. 

Velocity components. - Velocity components of the eddy-flow solution 
are shown in figures 6 to 8. These velocity components are directly 
related by equation (3) to the local partial derivatives of the velocity 
potential ^ given in figure 5. For the impeller geometry being inves- 
tigated, all these eddy-flow velocity components would be neglected by 
axial- symmetry- type solutions (reference l) . 
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Lines of constant (Q^) ^ on the meridional planes are shown in 

figure 6. This velocity component of the eddy-flow solution is tangent 
to the streamlines (constant q), and therefore to the velocities, of the 
axially symmetric flow through the stationary impeller. The velocity 
component (Qi)t has maximum values on the "blade surfaces (planes A 
and A') and is zero on the center plane E. This velocity component 
also becomes zero upstream and downstream of the impeller proper. 

Lines of constant velocity component (Qi)t| are shown on the 

meridional planes in figure 7. This velocity component of the eddy-flow 
solution is normal to the streamlines, and therefore to the velocities, 
of the flow through the stationary impeller. The velocity component 
(Q-^)^ has maximum values on the blade surfaces and is zero on the center 

plane E. This velocity component must also be zero along the hub and 
shroud boundaries, and becomes zero downstream of the impeller. Note 
that the velocity component (Qp)^ would be completely neglected in two- 

dimensional solutions on surfaces of revolution (reference 4) generated 
by streamlines of axial- symmetry-type solutions. 

Lines of constant tangential velocity component (W^)g relative to 

the impeller are shown in figure 8. For the impeller geometry investiga- 
ted, this velocity component has maximum values on the center plane E 
and is zero on the blade surfaces. Negative values of (Wp)© indicate 

flow across the meridional planes in the direction opposed to impeller 
rotation (into the page), and positive values of (Wp)^ indicate flow 

across the meridional planes in the direction of rotation (out of the 
page). From continuity considerations the integrated weight flow into the 
page (exclusive of the fluid that remains in the diffuser) must equal the 
integrated weight flow out of the page. These integrations have been 
carried out for the center plane E and weight flows agree within 

2j|j- percent. This agreement indicates approximately the same accuracy as 

that obtained from the integrated weight flows across the center line in 
figure 2(b). Thus it seems reasonable to conclude that the error through- 
out the domain of the three-dimensional solution is not greater than that 
indicated by the velocities in figure 2(b). 

Path lines . - Path lines of fluid particles on the passage surfaces 
are shown for the eddy-flow solution in figure 9. The fluid remains in 
the impeller passage and rotates in the opposite direction to that of the 
impeller . 

Solution for Flow Through Stationary Impeller Passage 

Flow through a stationary impeller with straight blades has zero 
tangential velocity and is equivalent to flow through the annulus formed 
by the hub and shroud surfaces. In reference 1, it is shown that for 
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incompressible flow the distributions of stream function and velocity 
components in the meridional plane for flow through an annulus are the 
same as the distributions for axially symmetric flow through a rotating 
impeller with an infinite number of straight impeller blades . . Therefore, 
example II of reference 1 is used in this report as the solution for flow 
through the stationary impeller passage. 

Streamlines. - Streamlines for flow through the stationary impeller 
passage are shown in figure 10. These lines are also the n coordinates 
(fig. 4(b)) used in the relaxation solution for the eddy flow. 

Velocity distribution. - For flow through the stationary impeller 
passage, (Q 2 ) 0 , (W 2 ) 0 , and (Qg)^ are zero. Lines of constant velocity 

(Qg) | (equal to Q 2 ) are shown on a meridional plane in figure 11. As 

for the eddy-flow solution, this velocity is expressed as a ratio of the 
tip speed of the rotating impeller, and the solution presented was 
obtained for (Q 2 )| equal to 0.3429 far upstream of the impeller proper. 

The distribution of (Q 2 )| is same ^ or a H meridional planes. 

Flow direction . - Lines of constant flow direction a 2 are shown on 
a meridional plane in figure 12. These values of a 2 can be used to com- 
pute the velocity components Qp and equation (10). 

Solutions for Flow Through Rotating Impeller Passage 

Solutions for various ratios of flow rate to impeller tip speed are 
obtained by superposition of various percentages (k in equation (5)) of 
the through-flow solution on the eddy-flow solution. Either the veloc- 
ity potential or the velocity components may be superposed. Two solu- 
tions are presented for flow through the rotating impeller with values 
of the axial inlet velocity Qz (equal to (Q 2 )| ) upstream of the 

impeller equal to 0.1372 and 0.3429, that is, for k equal to 0.4 and 
1.0, respectively. 

Solution for k = 0.4. - Path lines of fluid particles on the sur- 
face s — 3f"Th^ - impell^r - channel are shown in figure 13 for 40 percent of 
the through-flow solution superposed on the eddy-flow solution. Path 
lines on the hub and on the blade surface faced in the direction of rota- 
tion are shown in figure 13(a); path lines on the shroud and on the blade 
surface opposed to the direction of rotation are shown in figure 13(b). 

A composite plot of these path lines is shown in figure 13(c). 

For this solution the flow rate through the rotating impeller is not 
sufficient to eliminate (by superposition) all the reverse flow result- 
ing from the negative velocities (Qq)| the eddy-flow solution (see 

fig. 7(a))* This condition corresponds to the eddy flow that is attached 
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to the face of the hlade in the direction of rotation for two-dimensional 
solutions on surfaces of revolution (references 4 and 6, for example). 
Unlike the two-dimensional solutions, however, the fluid in the reverse 
flow of the three-dimensional solution does not remain in the impeller 
hut eventually leaves as indicated hy the spiral path lines emanating 
from the stagnation point on the huh of the impeller. 


The locus of stagnation points indicated on the hlade surface in 
figure 13(a) corresponds to the downstream stagnation point associated 
with the eddy flow of a two-dimensional solution. For the three- 
dimensional solution in figure 13(a), upstream stagnation points occur 
at the huh and shroud only. However, along the dot -dash line between 
these stagnation points the velocity component Q ^ is zero so that this 

line corresponds to the upstream stagnation point associated with the 
eddy flow of a two-dimensional solution. Path lines on the shroud sur- 
face in figure 13(h) converge to the upstream stagnation point. This 
convergence indicates that, as the path lines approach the stagnation 
point, the fluid leaves the shroud surface and passes into the interior 
of the passage. 


Solution for k = 1.0. - Path lines of fluid particles on the sur- 
faces of the impeller channel are shown in figure 14 for 100 percent of 
the through-flow solution superposed on the eddy-flow solution. The 
conditions for this solution are the same as those for the two-dimensional 
solutions given in references 1 and 4. Path lines on the huh and on the 
■blade surface faced in the direction of rotation are shown in figure 14(a) 
path lines on the shroud and on the hlade surface opposed to the direction 
of rotation are shown in figure 14(h). A composite plot of these path 
lines is shown in figure 14(c). 


COMPARISON OF TWO- AND THREE-DIMENSIONAL SOLUTIONS 

The results of the three-dimensional solution are compared with two- 
dimensional solutions on the mean passage surface extending from huh to 
shroud, on the mean surface of revolution, and on the shroud surface. 

O nl y the eddy-flow solutions are compared because the contribution of 
through flow to the velocity components is the same for both the two- 
and three-dimensional solutions. Thus the velocity components to he 
compared are components of the perturbation velocity caused hy the rota- 
tion of the impeller, and the relative importance of errors in these 
components is reduced when the known, primary through flow is added. 

Mean passage surface. - Because, for the type of impeller geometry 
investigated, the three-dimensional eddy flow has no velocity components 
in the mean passage surface extending from huh to shroud (plane E, 
fig. 4(a)), the velocity components in this plane are solely determined 
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"by, and therefore agree with, the axial- symmetry two-dimensional solu- 
tion. The axial- symmetry solution, however, completely neglects the 
relative tangential velocity W g, which for the three-dimensional solu- 
tion has maximum values on the mean plane (fig. 8(d)). 

For impellers with curved blades, the relative tangential velocity 
component reaches maximum values on a mean flow surface between the 
blades. If, as for high-solidity blade rows, this surface is not much 
different from the geometric mean surface between blades, then the flow 
is nearly two dimensional on the mean passage surface and the flow on 
this mean surface is approximately described by axial -symmetry solutions 
like those of reference 1. This conclusion is reached by Ruden in 
reference 8. 

Mean surface of revolution . - The velocity components Q| and Wg 

for the two- and three-dimensional solutions are compared on the mean 
surface of revolution in figures 15 and 16. The agreement for Q| in 

figures 15(a) and 15(b) is excellent, and the agreement for Wg in fig- 
ures 16(a) and 16(b) is also excellent near the impeller tip, although 
the two-dimensional solution (fig. 16(b)) introduces relatively small 
positive values of Wg not found for the three-dimensional solution in 
the region upstream of the contour line for Wg equal to zero. 

The slip factor, defined as the ratio of average absolute tangential 
velocity at the impeller tip to the tip speed of the impeller, depends on 
the distribution of Wq at the impeller tip and is equal to 0.7892 for 
the three-dimensional solution compared with 0.8142 for the two- 
dimensional solution on the mean surface of revolution (reference 4). 

The velocity component of the three-dimensional solution is 

plotted in figure 17. This velocity component is normal to the mean sur- 
face of revolution and is completely neglected by the two-dimensional 
solution. 

Shroud. - The velocity components Q| and Wg for the two- and 

three-dimensional Solutions are compared on the shroud surface in fig- 
ures 18 and 19. (The two-dimensional solution on the shroud surface was 
obtained from correlation equations, developed in reference 4, using, for 
"standard values" of velocity, the velocities of the two-dimensional, 
eddy-flow solution on the mean plane.) The agreement for in fig- 

ures 18(a) and 18(b) is good, but the agreement for Wg in figures 19(a) 

and 19(b) is poor, except in a limited region near the impeller tip. 

From simple physical considerations the agreement for and Wg on 

the hub is expected to be similar to the agreement on the shroud, except 
that Wg will have large positive values instead of the large negative 
values on the shroud. On the hub and shroud surfaces the velocity com- 
ponent Q-q is zero for both the two- and three-dimensional solutions. 
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Summary of comparisons . - A summary of the comparisons between the 
two- and three-dimensional - solutions that are discussed in this report 
is given in the following table: 


Velocity 

component 

Agreement 

Mean passage 
surface 

Mean surface 
of revolution 

Hub or shroud 
surface 

Q e 

good 

good 

good 


poor 

good 

poor 


good 

poor 

good 


It is concluded that on the flow surfaces investigated the velocity com- 
ponents W q , and agree for the two- and three-dimensional solu- 
tions discussed in this report, except: (l) W q on the hub, shroud, and 

mean passage surface, and (2) on the mean surface of revolution. 

If quasi -three -dimensional solutions are obtained by the proper 
combination of two-dimensional axial-symmetry and blade-to-blade solu- 
tions (reference 5), good agreement with the exact three-dimensional 
solution is indicated by good agreement on all surfaces of revolution. 
This agreement has already been discussed for the hub, shroud, and mean 
surfaces of revolution. For intermediate surfaces the table of compari- 
sons suggests that the agreement will always be good for Q|, will be 

progressively better for W q as the mean surface of revolution is 
approached, and will be progressively better for as the hub and 

shroud are approached. Because is the velocity component of prime 

importance, it is concluded that, at least for the type of impeller 
investigated in this report, two-dimensional solutions can be combined 
to describe the three-dimensional flow in rotating impeller passages with 
sufficient accuracy for engineering analyses. 


SUMMARY OF RESULTS AND CONCLUSIONS 

A solution is presented for three-dimensional, incompressible, non- 
viscous, potential flow in a rotating impeller passage with zero through 
flow. The solution is obtained for a conventional impeller with straight 
blades but with the inducer vanes removed and the impeller blades 
extended upstream parallel to the axis of the impeller. By superposition 
of solutions two additional examples are obtained for different flow 
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rates through the rotating impeller. Of particular interest is the fact 
that at low compressor flow rates the fluid in the reverse or eddy-flow 
region does not remain permanently in the impeller passage, as is the 
case for two-dimensional solutions on surfaces of revolution, hut, after 
spiraling around, eventually leaves the impeller . In other respects the 
three-dimensional solutions are compared with corresponding two- 
dimensional solutions and it is concluded that, at least for the type 
of impeller geometry investigated, two-dimensional solutions can he 
combined to describe the three-dimensional flow in rotating impellers . 
with sufficient accuracy for engineering analyses. In particular it is 
concluded that: 

1. On the mean surface of revolution the velocity components, 
except the component normal to the surface, agree for the two- and 
three-dimensional solutions. 

2. On the huh and shroud surfaces the relative tangential velocity 
component does not agree for two- and three-dimensional solutions, hut 
the other velocity components do. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, July 3, 1952 
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APPENDIX A 
SYMBOLS 

The following symbols are used in this report: 


A, A' ,B,B' , . . . 

. meridional planes (fig. 4(a)) 

k 

percentage of through-flow solution 

Q 

absolute velocity, expressed as ratio of impeller tip 
speed 

R, 0,Z 

cylindrical coordinates (fig. 1(a)), linear coordinates 
expressed as ratios of impeller tip radius 

W 

relative velocity, expressed as ratio of impeller tip 
speed 

02 

angle, figure l(b) and equation (8) 


velocity potential and stream function, respectively, for 
incompressible flow through hub-shroud annulus; used as 
coordinate system in £q-plane, equations (Bl) and (B2) 

cp 

velocity potential, equation (3) 


Subscripts: 


1 

rotating impeller with zero net through flow (eddy-flow 
solution) 

2 

stationary impeller with through flow (for numerical 
example, through flow is such that Qz equals 0.3429 
upstream of impeller) 

R,9,Z, 

components in R,0,Z, £,t] directions, respectively 
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APPENDIX B 


TRANSFORMATION FROM RZ- TO |t) -PLANE 

It is convenient for purposes of solution by relaxation methods to 
transform the RZ-plane to one on which the coordinates are the stream 
function T) and the velocity potential | for flow through the com- 
pressor annulus without blades. The stream function t\ satisfies the 
continuity condition if defined as 


5r = R(Q 2 )z 

Sz = " R(Q2)r 


(Bl) 


and the velocity potential | satisfies the irrotationality condition 
if defined as 



(B2) 


In terms of the transformed r ) 

tives of equation (4b) become 


coordinates, the partial deriva- 

\ 


dcp dcp dcp dr) 
Sr ~ Sf Sr + 5ti Sr 


S 2 cp _ Sfcp f S i ] 2 

3e 2 * a { 2 w 


+ 2 


S 2 cp S| St] 
S|Sq Sr Sr + 




(B3) 


Sep Sep S£ Sep Sti 

Sz = S| Sz + S^ Sz 


Sz 2 s $ 2 



S 2 cp SE_ Sn 

+ 2 S^ Sz Sz + 
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From equations (Bl) through (B3) , equation (4b) becomes 

/ 2 r\ 2\ ( d 2 cp -d 2 d 2 cp\ 1 d 2 qp % ^z\ dcp 

(<iR + 152 )2 (ajf + R i $) + ? 5? + (55- + t + ZTJz sf 


2 (q 2 )z + R (sr " zr ) 2 


S?-o 


(B4) 


But, 


(Or 2 + %, 2 ) 2 ~ Q 2 2 


and. from continuity 


dQ R Q R SQ z \ 

5r~ + T + 3z"/2 = 0 


and for irrotational flow 

/dQ z ^Q R \ 

V5iT - 5 z"/2 = 0 

so that equation (B4) becomes 


Q 2 ‘ 


dfcp 

H 2 


R 2 Q 2 2 


S 2 cp _1_ c^cp 
dr) 2 + R 2 d0 2 


2 (^2 )z 5^ = 0 


( 6 ) 


Equation (6) is the partial differential equation for the distribution 
of cp in the £, t \>9 coordinate system. 

This new coordinate system introduces two new velocity components 
and which are related to the radial and axial velocity com- 

ponents by equation (9). Combining equations (3), (9), and (Bl) to 
(B3) gives 


Q* = Q 2 §j (7a) 

S = rq 2 5? 


(7b) 
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(a) Cylindrical coordinates. 

Figure l. - Coordinate systems relative to impeller, and absolute velocity 
components . All quantities are dimensionless. Linear coordinates are 
measured in units of impeller tip radius; velocity components are 
measured in units of impeller tip speed. 
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(b) £,ti coordinate system in RZ-plane. 

Figure 1. - Concluded. Coordinate systems relative to impeller, and absolute 
velocity components. All quantities are dimensionless. Linear coordinates 
are measured in units of impeller tip radius; velocity components are 
measured in units of impeller tip speed. 
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Direction of rotation 



(a) Comparison of graphically constructed path lines (three-dimensional solution) with streamlines 
obtained from stream function of two-dimensional solution. 


Figure 2. - Comparison of results obtained from two- and three-dimensional solutions. Plane normal 
to impeller axis in region of uniform axial velocity far upstream. 
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Direction of rotation 


+ W s 



(b) Lines of constant velocity relative to rotating impeller. 

Figure 2. - Concluded. Comparison of results obtained from two- and three-dimensional solutions. 
Plane normal to impeller axis in region of uniform axial velocity far upstream. 
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Figure 3. - Hub- shroud dimensions of impeller for numerical 
examples. Vaneless diffuser; straight impeller blades 
extended far upstream parallel with axis of impeller. 
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Impeller axis 
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Impeller axis 
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Impeller axis 
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(b) Merldlortal plane B. Values of (Ql)$ are same on plane 
of opposite sign. 

Figure 6. - Continued. Lines of constant (Qi)$, velocity component of eddy-flow 
solution tangent to velocity of through-flow solution. 
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Impeller axis 
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Impeller axis 
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Impeller axis 
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(c) Meridional plane 


Values of are same on plane C* but 

of opposite sign. 

Lines of constant (Q^)^ velocity component of eddy-flow 


Figure 7. - Continued. 

solution normal to velocity of through-flow solution and lying In meridional 
plane . 
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Impeller axis 
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Impeller axis 
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Impeller axis 
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(d) Meridional plane E. 

Plgiire 8. - Concluded. Lines of constant (W^g, 

of eddy-flow solution, relative to rotating Impeller, 
are directed Into page. 



Impeller axis 
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Figure 9. - Path lines or riuld particles on surfaces of 
rotating Impeller channel with zero net through flow. 
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(b) Blade surface opposed to direction of 
rotation, and shroud. 

Figure 13. - Continued. Paths of fluid particles on 

surfaces of rotating channel. Forty percent of through- 
flow solution superposed on eddy-flow solution (k - 0.4); 
$2* 0.1372 far upstream. 
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(c) Composite plot. 

Plgurc* 13. - Concluded. Paths of fluid particles on sur- 
faces of rotating channel. Forty percent of through-flow 
solution superposed on eddy-flow solution (k - 0.4)j Q z , 
0.1572 far upstream. 
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(c) Composite plot. 

Plgure 14. - Concluded, i’aths of fluid particles or. sur- 
faces of rotating channel. One hundred percent of through 
flow solution superposed on eddy-flow solution ,'k - 
Qy, 0.3429 far upstream. 
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Figure 15. - Mean surface of revolution showing velocity 
component in direction of through flow. Eddy- 

flow solutions. 
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Figure 16. - Mean surface of revolution showing tangential 
velocity component W@ relative to rotating Impeller. 
Eddy-flow solutions. 
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